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ABSTRACT 

Observations show that star formation in galaxies is closely correlated with the abun- 
dance of molecular hydrogen. Modeling this empirical relation from first principles 
proves challenging, however, and many question regarding its properties remain open. 
For instance, the exact functional form of the relation is still debated and it is also 
unknown whether it applies at z > 4, where CO observations are sparse. Here, we 
analyze how the shape of the star formation - gas relation affects the cosmic star for- 
mation history and global galaxy properties using an analytic model that follows the 
average evolution of galaxies in dark matter halos across cosmic time. We show that a 
linear relation with an H2 depletion time of ~ 2.5 Gyr, as found in studies of nearby 
galaxies, results in good agreement with current observations of galaxies at both low 
and high redshift. These observations include the evolution of the cosmic star forma- 
tion rate density, the z ~ 4 — 9 UV luminosity function, the evolution of the mass - 
metallicity relation, the relation between stellar and halo mass, and the gas-to-stellar 
mass ratios of galaxies. In contrast, the short depletion times that result from adopting 
a highly super-linear star formation - gas relation lead to large star formation rates, 
substantial metal enrichment (~ 0.1 Zq), and low gas-to-stellar mass ratios already 
at z > 10, in disagreement with observations. These results can be understood in 
terms of an equilibrium picture of galaxy evolution in which gas inflows, outflows, and 
star formation drive the metallicities and gas fractions toward equilibrium values that 
are determined by the ratio of the gas accretion time to the gas depletion time. In 
this picture, the cosmic modulation of the accretion rate is the primary process that 
drives the evolution of stellar masses, gas masses, and metallicities of galaxies from 
high redshift until today. 

Key words: galaxies: evolution - galaxies: ISM - galaxies: high-redshift - galaxies: 
star formation 



1 INTRODUCTION 

The relation between the star formation rate (SFR) per unit 
area, Esfr, and the gas surface density (of either molecular 
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I Rahman et al.l l201ll: iFeldmann et al.l l201ll: iRahman et al.l 
l2012l : lLada et al.ll2012l : IFeldmann et alj|2012h . This relation 
holds over several orders of magnitude in Esfr with a 
relatively small scatter if measured on ~ kpc sized patches 
of galaxies or for whole galaxies. Yet, despite these efforts, 
we still do not quite understand how this correlation comes 
about and, embarrassingly, what the precise functional 
form of this relation is. 

The reason for the latter is that neither SFRs nor molec- 
ular hydrogen (H2) masses are directly accessible in extra- 
galactic observations. SFRs are inferred from a combination 
of ultraviolet (UV) and infrared (IR) observations, in order 
to capture the energy output of both obscured and unob- 
scured star formation. Cold molecular gas is measured indi- 
rectly from emission lines of tracer molecules, such as car- 
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bon monoxide (CO). While observational studies of nearby 
galaxies and galaxies out to z ~ 3 seem to show that the re- 
lation between Ssfr and E h, is close to linear (jBigiel et al.l 
120081 ; iDaddi et al1l2010al lbl; iGenzel et alll2010l ), there is no 
clear theoretical understanding why this is and, furthermore, 
there are potential concerns that systematics in the observ- 
ables bias these results. For instance, it has been suggested 
that the star formation - gas relation is super-linear with a 
power- law exponent of ~ 1.4 on kpc scales, but that IR emis- 
sion unrel ated to rec e nt sta r formati on reduces the i nferre d 
slope, see iLiu et al.l i|201ll ). but cf. iRahman et~aH (|201ll ). 
Alternatively, systematic variation of the CO/H2 conver- 
sion factor with CO intensity, SFRs, or gas surface den- 
sity could potentially disgui se a highly super-linear relatio n 
an d make it appear linea r l|Naravanan et all 1201 ll . l2012bl ) , 
cf. lFeldmann et al.l (|2012l ). Finally, the fitting itself has been 
targeted as a potential s ource of bias. A Bay esian re-analysis 
of the data presented bv lBigiel et al.l (|2008T > results in tenta- 
tive evidence that the star formation - gas relation is mildly 
sitb-linear a nd shows substanti al variations from one galaxy 
to the next (|Shettv et al.ll2012f ). 

Theoretical explanations of the star formation - 
gas connectio n predict a variety of non-linear relation- 



ships, s e e e.g., Krumholz fc McKee (120051); Krumholz et alj 



( 2009bh: lOstrilcer et al 



I 20ld ); lOstriker fc Shettvl J201ll) 



For instance, Krumholz et al.l 
the power-law exponent is mildly 
1. 33) at high gas su r face densities (> 100 



IKrumholz et al.l <|2012T ). 
(20093) estimate that 
non-linear (" 

Mq pc~ 2 ), while lOstriker fc Shettvl (120111 ) predict a much 
steeper, quadratic relation at such gas surface densities. This 
gas surface density regime is clearly highly relevant because 
it corresponds to typical gas surface densities of star fo rm- 
ing galaxies at z > 1 and higher (e.g. jGenzel et al"1l2010] . see 
also Section [3]2]). A significantly super- linear star formation 
- gas relation has also been suggested by recent work based 
on (non-cosmological) galaxy evolution simulations that in- 
clude a ~ pc-scale modeling of v arious stellar feedback chan- 
nels l|Hopkins et al.ll201ll . 12012?) . 

In this paper we take a very different approach to 
put constrains on the functional form of the star forma- 
tion - gas relation. We study how global galaxy proper- 
ties vary depending on the chosen star formation - gas 
relation and compare our model predictions with obser- 
vations. In particular, we show that the evolution of the 



cosm i c SFR density (SFR D; iLillv et al.ll"l996l; iMadau et al 
' Steidel et ail Il999l) at z ~ 4 - 10 (iBouwens et al 



1996 



2007 



2012 j ; iBehroozi et al.l l2012bl ) is highly sensitive to 
the adopted form of the star formation - gas relation, and, 
hence, that it can be used to test its functional form. 

The paper is organized as follows. In Section [2] we de- 
scribe the analytic model which underlies our predictions 
and present the star formation - gas relations that enter our 
analysis. In Section [3.11 we demonstrate that the model re- 
produces the observed stellar mass - metallicity relation, the 
stellar mass - halo mass relation, and the gas-to-stellar mass 
ratios of galaxies reasonably well if a linear star formation - 
gas relation is chosen. We also show that predictions based 
on a quadratic star formation - gas relation are in conflict 
with observations. Next, in Section |3. 21 we present how the 
cosmic star formation history is affected by the functional 



form of the star formation - gas relation. We discuss the 
origin of this result in Section[3]3] In Section [3]4] we put our 
findings in a broader context, showing that galaxy evolution 
can be thought of as an evolution along a line of fixed points 
(equilibria) driven by the cosmic modulation of the gas ac- 
cretion rate. Potential caveats of our study are presented in 
3.51 We summarize our findings in Section [4] and close the 
paper with a comparison to previous work. 



2 THE MODEL 

The analytic model that is used in the remain der of the pa- 
per fo llows closely the approach pres ented bv iBouche et al.l 
(|20ld) and lKrumholz fc Dekel i|2012l ). At its most basic level 
the model tracks the average evolutiorQ of the dark matter, 
gas, and stellar content of halos using theoretically or obser- 
vationally motivated accretion and transformation rates for 
the different mass components. For the convenience of the 
reader we begin by summarizing the main equations that 
enter the model and refer the reader to IKrumholz fc Dekell 
|2012l ) for more details. We then discuss the various star for- 
mation models that we explore in this work and conclude the 
section with a short descriptions of the actual set-up of the 
model. In the following, we adopt co smological parameter s 
consistent with WMAP 7-year data (|Komatsu et al.ll20lih : 
n m = 0.27, n b = 0.046, n A = 0.73, h = 0.7, and a s = 0.81. 



2.1 Basic equations 

Accretion rates: The total (baryon plus dark matter) ac- 
cretion rat e at time t o nto a halo of a given mass 
Afhaio is (jNeistein et al ] 120061 : iNeistein fc Dekell [2008; 
IKrumholz fc DekellbOll ) 



M ha i = 1.06 x 10 i2 A/« 







M h 



1O 12 M 



D(t) 
D(ty 



(1) 



D(t) is the linear growth factor normalized to unity at the 
present epoch. The (average) evolution of the mass of a halo 
can be computed from its initial mass and the time- and 
mass-dependent accretion rate (Q]). The corresponding ac- 
cretion rates of the gas and stellar components are 



M gi i n — £in/g,in/b -Mhalo, 
M,,in = (1 - /g,in)/b Afhalc 



(2) 
(3) 



Here, ei n specifies which fraction of the accreted gas reaches 
the galaxy at the center of the halo, /b = fib/^m is 
the cosmic baryon fraction, and / g ,i n is the gas fraction 
(in units of /b) of the accreted matter. We assume that 
accretion of gas onto the disk is completely quenched 
(ei n = 0) when the halo mass is above the virial shock 
mass M s ho ck (a)/10 12 M B = max(2,10 12 M B /(3M t ,Ps)) 
l|Dekel et al.ll2009l ; IKrumholz fc Dekel2012l ). Here, M», PS (z) 
is the Press-Schechter mass at the given redshift, see 



1 Since the model does not use merger trees, it only follows the 
mean evolution of halos. For the quantities studied in this pa- 
per (the cosmic star formation history and the gas fraction of 
galaxies), such an average description is adequate. 
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iMo fc White] l|2002l ). Accretion is also quenched when the 
halo mass is so low that photo-ionization limits the baryon 
fraction in the halos to 50% of /b or less (|Hoeft et al.ll200rj| ; 
lOkamoto etHl l200Sfi . Specifically we assume ei n = if 
the halo mass i s belo w the mass M c (z) given in Fig. 3 of 
lOkamoto et~ai1 (|2008T ). We adopt M c = for z > 9. We 
assume that accretion is not quenched (ei n = 1) if the halo 
mass lies between M c and M B h ock. W e compute / g ,i n as de- 
scribed bv lKrumholz fc Dek cl (2012). In brief, we first com- 
pute which fraction of the accretion rate is due to accreted 
halos of a given mass. The model provides us directly with 
the gas fraction of a halo of a given mass and we then cal- 
culate / gi i n as the accretion-weighted mean gas fraction (in 
units of /b) of all accreted halos plus an additional non-halo 
contribution of 20%. 

Evolution of the gas and stellar mass: Besides accretion, 
the baryonic components are subject to additional physi- 
cal processes, such as star formation, stellar mass loss, and 
galactic outflows which we consider next. Specifically, we 
assume that the gas and stellar masses change according to 



M e ,in- (l--R+e out )M» 



,form j 



M, = M», in + (1 - fl)M*, for: 



(4) 
(5) 



M*,f rm is the SFR of the galaxy (see below). R is the frac- 
tion of the stellar mass of a single stellar population that is 
recycled to the ISM vi a, e.g., stellar w inds. For the IMF 
assumed in this work (|Chabrierl 120051 ) . R = 0.46 is ap- 
propriate. The model makes the simplifying assumption of 
instan t recycling which is , how ever, adequate for our pur- 
poses l|Krumholz fc Dekell liiilh . The parameter e out is in- 
cluded to model l a rge s cale galactic outflows. We follow 
iKrumholz fc Dekell (|2012T l and assume e ut — 1, i.e., a gas 
outflow ra te equal to the SFR, based on the estimates by 
l|Erbll2008l ). We discuss the impact of varying this and other 
model parameters in Section |3, 5 1 

The gas disk: For the modeling of the star formation 
(see below) we need to know how the gas is distributed in 
the galaxy. We make t he ansatz that the ga s forms an ex- 
ponential disk, see e.g., iBigiel fc Blit3 (|2012h . 



E g (r) = E g (0)exp-^. 



(6) 



We assume that the scale radius i?d is se t by the initia l, 
pre-collapse angular momentum of the gas (|Mo et al.ll 19981 ) . 
Specifically, 



i? d = 0.5Ai? vir , (7) 

where i? vi r(M halo , z) and A = 0.07 ([Krumholz fc Dekel20l"3 ) 
are the virial radius and (median) spin parameter of the 
halo, respectively. The prefactor takes into account that the 
half- mass radius of an exponential disk is ~ 1.7J?d- The pre- 
dictions of equation Q are in re asonable agreement with 

l|2012l ) find a mean scale radius 



cqi 

observations. IBigiel fc Blitz 
of 6.1±0.7 kpc (with significant galaxy-to-galaxy variations) 



for a sample of roughly MW-like galaxies in the local Uni- 
verse, while equation predicts Rd ~ 9 kpc for a 10 12 
Mq halo at z = 0. Reducing A to a more traditional value 
A ~ 0.04 - 0.05 dBarnes fc Efstathioulll987l ; iBullock et all 
l200ll ; iDutton et al. 2011 ) , would lead to even better agree- 
ment, see also Section |3. 5 1 



Molecular hydrogen: The abundance of molecular gas 
plays an important role in many of the star formation models 
that we consider in the paper (see below). We compute the 
H2 surface d ensity £h 2 = /h 2 E g following IKrumholz et al.l 
(2008. l2009ah : 



/ H2 (E g , Z) = 



if s < 2 



0, otherwise 

ln(l + 0.6x + 0.01x 2 ) 
0.6r c 

^(1 + 3.1(Z/Z e f 365 ) 
E e 



3T 
4.1 

320 c 



gem- 



;Z/Z Q , c =5. 



(8) 

(9) 
(10) 
(11) 



Here, Z — Mz/M s is the average metallicity of the gas and 
we effectively assume that there are no strong metallicity 
gradients. We note that gas disks in nearby galaxies often 
show a moderate radial metallicity gradient of the order of ~ 
0.05 — 0.1 dex per kpc (|Zaritskv et al.| [l994'). corresponding 
to a factor ~ 1.5 — 4 decrease in oxygen abundance from 
the c enter to the optical edge of the disk (jMoustakas et all 
2010). We discuss how metallicity gradients affect our results 
in Section 13,51 

Metallicity evolution: The dependence of the H2 abun- 
dance on the metallicitjjfl of the ISM requires that we keep 
track of the build-up of metals in each galaxy. In the in- 
stantaneous recycling approximation and ignoring inflows 
and outflows the total metal production rate is propor- 
tional to both the yield y of the stellar population and the 
star formation rate, i.e., Mz.tot = 3/(1 — R)M, t torm- Met- 
als that are locked in stellar remnants do not contribute 
to the metallicity of the ISM, which therefore evolves ac- 
cording to Mz = y(l - fl)M«,fo r m - Z (l - i?)M*, form . As 
suggested by IKrumholz fc Dekell l|2012l ) we add three pro- 
cesses to the basic picture. First, a fraction C(^haio) = 
0.9 exp(— Mhaio/3 x 1O 11 M0) of the newly produced metals 
(predominantly in a hot, shocked phase) is removed from 
small halos by supernovae driven outflows. Second, cold and 
warm gas of the ISM and the associated metals are lost by 
galactic winds at a rate e ut(Mz/M g )Af, i f orm . Third, inflows 
from the IGM add metals at the rate ZiGM-^g.in- Hence, the 
total metal mass of a galaxy changes according to 

M Z = [2/(1 - R)(l - ~ (1 - R + £out)M Z /Mg] inform 

+ ZiqmM sM . (12) 

A single pair instability supernova ca n enrich halos at high 
redshifts to metallicities ~ 2 x 10" 5 dWise et al.ll2012h . We 
therefore adopt the value Z\gm = 2 x 10 -5 as the metallicity 
of both the initial gas and the gas accr e ted on to halos. 



Differences to Krumholz & Dekel \20 lA) : Our imple- 



mentation of_tfT«_jriodel_fiollows_closely the approach pre- 
sented by IKrumholz fc Dekell (|20 12T ) . The only noteworthy 
differences are the following. 



2 More precisely, it is the dust-to-gas ratio that regulates the 
abundance of H2. We assume that the dust-to-gas ratio scales 
proportion al with metal l icity a s observed for galaxies enriched to 
Z > 0.1 Z Q iLerov et al] (l201lh . 
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• We added a mass threshold below which gas accretion 
is quenched due to photo-ionization. 

• The metallicity evolution of the ISM (Equation I12|) 
accounts for the locking-up of metals in stellar remnants. 
Math ematically this amounts to replacing e out in Equation 
23 of iKrumholz fc Dekell (|2012l ) with e ou t + (1 — i?) and, 
hence, is equivalent to a change of the parameter e ou t. We 
show in Section 13.51 that varying e ou t by a factor of a few 
does not strongly affect the SFRDs at z > 2. 

• We use a variety of different star formation - gas rela- 
tions (see next section). 



2.2 Star formation 

We calculate the SFR of a galaxy by converting a gas surface 
density (either the total gas or the H2 surface density) into 
a SFR surface density with the help of an empirically moti- 
vated or theoretically derived star formation - gas relation. 
The SFR of a galaxy is then simply the integral of the SFR 
surface density over the area of the gas disk of the galaxy. 

A variety of star formation - gas relations have been 
suggested based on theoretical modeling and/or observa- 
tions. The relationship often takes the form of a power-law 
with an exponent between approximately 1 (linear) and 2 
(quadratic). In this paper we primarily focus on the two 
limiting cases out of the many possibilities for the following 
reasons. First, comparing a linear with a quadratic relation 
makes the qualitative changes in global galaxy properties 
more apparent that arise from the functional form of the 
star formation - gas relation. Second, we find that rela- 
tions with exponents intermediate between the linear and 
the quadratic case result in corresponding intermediate pre- 
dictions. Hence, much can be learned already from a study 
of the more extreme cases. Finally, both the linear and the 
quadratic star formation - gas relation are supported by 
various lines of evidence. It is thus clearly of high interest 
to test their specific implications for the evolution of galaxy 
properties. 

An approximately linear relation is empirically moti- 
vated by observationa l studies of the S s fr — Eh 2 relation 
in the local universe (Bigicl ct al. 2008, 1201 ll ) and out to, 
at least, z ~ 2 |Genzel et al. 2010T). The depletion time 
in nearby galaxies is ~ 2 — 3 Gvr (|r3igiel et all 120111 ). but 
whether this is also true at high redshift is not so clear, 
see, e.g., iTacconi et al.l (|2012h . High redshift samples typi- 
cally include only galaxies with large star formation rates 
and gas masses, and, hence, in case of a slightly super-linear 
star formation - gas relation it is difficult to distinguish a 
change in the depletion time with redshift from a change 
in depletion time due to the increase in gas surface den- 
sity. Uncertainties in the CO/H2 conversion factor and the 
limited spatial resolution of CO observations further com- 
plicate an accurate measurement of the depletion time at 
high redshift. A quadratic relationship (at gas surface den- 
sities > 50 Mq pc -2 ) is expected based on analytic and 
numerical models of feed back-regulated star formation, e.g., 
lOstriker fc Shettvl ((2011). It has also been argued that a 
quadratic relationship is supported by current observations 
once the conversion between CO emission and H2 surface 



densit y is pr operly taken into acc ount (jNaravanan et ail 
l2012bl . but cf. lFeldmann et al.ll2012h . 

In addition to the mentioned linear and quadratic re- 
lationship we also explore how alternative star formation - 
gas relations would affect the cosmic star formation history. 
In the following E may refer to either the H2 surface density 
Eh 2 = /h 2 (Eg, Z)E g (the default), to the H2 surface density 
Eh 2 ,z = /h 2 (Eg, z?0)E g for gas of solar metallicity, or to 
the total gas surface density, E g , respectively. The full list 
of star formation - gas relations used in this paper is the 
following: 

• A linear relation w ith a constant depletion time of 2.5 
Gyr jBigiel et al.ll201lD , i.e., 



EsFR = 



(13) 



2.5 Gyr 

• A relation that is quadratic at gas surface densities 



^ 50 Mi 



2 P c 



jOstriker fc Shettvll201lh : 

O.lMoyr^kpc -2 



and linear at smaller surface densities 



, r w „, , I if Sg ^5OM pc- 

ESFR = < |_ M 0P C J 

[2^ ifE g <5OM pc- 

(14) 

The z > 3 predictions of this star formation - gas rela- 
tion are effectively indistinguishable from the predictions 
of a relation that extends the quadratic relationship to 
E g < 5OM0pc -2 , because surface densities of star form- 
ing galaxies at high redshift typically exceed 5OM0pc~ 2 , 
see Section [3.21 

• A relation with a moderately s ub-linear exponent 0.84 
as suggested bv lShettv et~ail (|2012T ) (with E = Eh 2 ) : 



10 M Q pc- 



084 10M Q pc- 2 



2.5 Gyr 



(15) 



• A moderately su per-linear relation suggested by 
IKrumholz et al.l (|2009bh (with E = Eh 2 ): 



EsFR = 



2.6Gyr 



85M Q pc- 
v 85Mqpc~ 



if Eg 
if E, 



< 85 M Q pc" 2 

> 85 M Q pc" 2 . 

(16) 

• The "unive r sal" s tar formation relation suggested by 
IKrumholz et al.l (|2012r i (with E = Eh 2 ): 



E S fr = 0.01 — , where tg = min(t ffiG MC, tg tT ), 
iff 



(17) 



and 



iff,GMC = 



rV4 



, Egmc = 85M©pc 



tff,T 



(18) 
(19) 



G'(E^ MC Eg) 1 / 4 

ZEi 
8 n 

In Equation (|18[) . a is the gas velocity dispersion on scales 
comparable to the disk scale height (few 100 pc). We set 
it to 8 km/s, adequate for disk galaxies in the local uni- 
verse ((Walter et alj [2008). Hi gh redshift galaxie s may have 
larger velocity dispersions (|Tacconi et al.|[2010h. but since 
typic ally tff,GMC > tff.T for these galaxies ( Krumholz et alj 
120121 ). the choice of a does not matter in this case. Equation 
(|19[1 includes the angular velocity of galactic rotation Q. We 
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Figure 1. Metallicity vs stellar mass. (Left) Model predictions based on a linear Esfr — Sh 2 relation, Equation (1 1 3 D - (Right) Model 
predictions based on the quadratic relation given by Equation l|14j l. Solid lines are the model predictions for z = 0.1, z = 2.2, z = 3.5, 
z = 6.0, and z = 9.0 (see legend). The dot-dashed line (m edian) and the shaded region (1-cr scatter) at the top show the observed 
low re dshift (z ~ 0.1) mass -metallicit y relation from SPSS l|Tremonti et alj|2004 ). The dotted line shows the results of a re-analysis 
of the iTremonti et al 1 d2004) study by iKewlev fc Ellison! J2008l l. Circles and squares ar e observational e stima tes for z ~ 2.2 (based on 
stacked spectra: lErb et al.ll2006}l and for z ~ 3.5 (indi vidual galaxies), respectively, by iMaiolino et all l|2008h . who use the same IMF 
and metallicity calibrations as iKewlev fc Ellison! [120081 1. Predictions based on a linear Ssfr — ^H 2 relation with a constant ~ 2.5 Gyr 
depletion time are in good agreement with observations of the mass-metallicity relation at z > 2. 



chose Q — Vn>t/(1.7 Rd) as our fiducial value. Here, 1.7 Rd 
is the half-mass radius of the gas disk and V ro t is the rota- 
tion velocity of the disk at 1.7 Rd- The latter we equate with 
V^nax,h, the maxim um circular velocity of a dark matter halo 
with NFW profile (|Navarro et al.lll997l ) that hosts the given 
galaxy. We determine Knax,h for a halo of given mass and 
at a given redshift using the fit s (conv erted to our definition 
of halo mass) by IPrada et al.l (|2012|) based o n high reso- 
lution N-body simulations ( Klypin et al.ll20l"ll ). The ansatz 
Vrot ~ Knax.h has a firm observational j ustification, at leas t 
for disk galaxies in the local universe (|Reves et al.l I2012T ). 
In compact star formation galaxies at high redshift, how- 
ever, it is conceivable that V r ot exceeds Knax,h- If the latter 
is true, we potentially overestimate tff^ and, hence, poten- 
tially underestimate the SFR at high redshift. As we will 
show in Section T3, 2 1 this only increases the conflict between 
this particular star formation - gas relation and observa- 
tions. 



2.3 Running the model 

We have implemented the model described in the previous 
sections as a GNU Octavjfl pro g ram. We adopt the fidu- 
cial choice of lKrumholz fc Dekell l|2012i ) (see their Table 1) 
for the model parameters, and note that most values are 



3 http:/ /www. gnu.org/software/octave 



justified based on either theoretical grounds (e.g., the recy- 
cling fraction R is a consequence of the adopted IMF) or on 
observations (e.g., /b, y). We discuss the robustness of our 
predictions to changes in the parameter values in Section 
1331 

We assume that each halo initially contains a fraction 
/b ~ 0.17 of its total mass in baryons. The large majority 
of the baryons is assumed to be in a gaseous component 
(90%) with an initial metallicity of Zjgm = 2 x 10~ 5 . The 
remaining 10% are locked up in a stellar component fo rmed 
during mergers at high redshift l|Rodighiero et al EuTlh . We 
checked that decreasing the initial stellar contribution from 
10% to 1% (and increasing the initial gas contributions cor- 
respondingly) does not change any of our results in a signif- 
icant way. We start the model at z — 30 and run it forward 
until z = in 300 steps spaced equally in values of the 
expansion factor a = 1/(1 + z). The model follows the con- 
current evolution of 300 halos of different mass ranging from 
3.2 x 10 5 M Q to 1.3 x 10 8 M G at z = 30 and from 1.4 x 10 7 
M Q to 8 x 10 15 M Q at z = in a logarithmic spacing. 



3 RESULTS 

In this section, we analyze the importance of the functional 
form of the star formation - gas relation for various galaxy 
properties and the cosmic star formation history in the con- 
text of the model of Section [2~Tl We will show that a strongly 
non-linear star formation - gas relation results in predic- 
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Figure 2. Gas-to-stcllar mass ratio in galaxies as function of their stellar mass. (Left) Model predictions based on a linear Esfr — 2h 2 
relation, Equation ( 1 1 3 1) - (Right) Model predictions based on the quadratic relation given by Equation 1141 . The predicted H2-to-stellar 
mass ratios are shown by solid lines, while the total gas- to-stellar mass ratio are shown by dotted lines. Different colors correspond to 
different rcdshifts as i ndicated in the legend . The circles (arrows) show > 3 — a detections (upper limits) of the H2 fraction of galaxies 
in the local universe jSaintonee et alj|201ll) . Circl es (2 ~ 0.5), square s (2 ~ 1.5), and stars (2 ~ 2.5) are measurements of cold gas 
masses for a sample of galaxies at higher redshift llMagdis et al. I l2012al) . These latter observations are based on dust infrared emission 
and are thus not affected by uncertainties related to the CO-H2 conversion. A linear Esfr — £h 2 relation predicts a large increase in 
the gas fractions up t o ~ 2 in good agreement with the observed evolution, followed by only a modest change toward higher redshift, cf. 
iMaedis et all l|2012bl ). 



tions that are in conflict with observations, while predictions 
based on a linear star formation - gas relation match the var- 
ious observational data better and in most cases remarkably 
well. 



if we had used a purely quadratic relation down to low gas 
surface densities, the mass - metallicity relation would not 
evolve between z = 9 and z — - in obvious disagreement 
with observations. 



3.1 Metallicities, stellar and gas masses 

In Fig. [1] we show the evolution of the stellar mass - 
metallicity relation as predicted by the model for a lin- 
ear and a quadratic Esfr — Eh 2 relation, see Equations 
(|13p and (|14|) with E = Eh 2 - We compa re these predic- 
tions with observations of gala xies at z ~ ( Tremonti et alj 
20041; iKewlev fc Ellisonl |200S| ). at z ~ 2 (lErb et al l 120061 ; 
Maiolino et alj|2008h . and at z ~ 3.5 (|Maiolino et al.ll200ct ). 
Predictions based on a linear relation are in excellent agree- 
ment with the observed mass - metallicity relation at 2 > 2. 
At low redshifts and low stellar masses the predicted metal- 
licities fall somewhat short of the measured ones, although 
we could reduce this discrepancy significantly (to less than 
0.2 dex) by decreasing the maximum value of C, (the expelled 
gas fraction see Section [2| from 0.9 to 0.8, which would also 
not strongly affect any of our other predictions. A quadratic 
relation strongly underpredicts the evolution of the mass - 
metallicity relation at high redshift. The evolution at z < 2 
is similar to the one for the linear star formation - gas rela- 
tion, because as those times typical gas surface densities in 
galaxies falls below ~ 50 M© pc -2 and the star formation - 
gas relation is effectively linear , see Equation (| 14[) . In fact, 



In Fig. [2] we compare the predicted gas-to-stellar mass 
ratios with observations spanning the redshift range z ~ 
— 3. A quadratic star formation - gas relation underpre- 
dicts both the total gas- and the molecular gas-to-stellar 
mass ratios. At z ~ 0.5 the model predictions are off by a 
factor ~ 3. The disagreement becomes even more severe at 
higher redshift. At z ~ 2 the observed gas-to-stellar mass ra- 
tios exceed the predicted ones by an order of magnitude. In 
contrast, we find good agreement between model predictions 
and observations for both the normalization and the redshift 
evolution of the gas-to-stellar mass ratio in case of a linear 
star formation - gas relation. For instance, the predicted 
fraction Mh 2 /(Mh 2 + M») is 39% at 2 = 1.5 and 54% at 
2 = 2.5, co mpared with 34% at 2 ~ 1.2 and 44% at 2 ~ 2.3 
as found bv lTacconi et all (|2010h based on CO observations 
of a sample of massive star forming galaxies. However, it 
has been argued that the observed gas masses are poten- 
tially overestimated, because of systematic variations of the 
CO/H2 conversi on factor with CO surfac e brightness that is 
unaccounted for lNaravanan et al.l (|2012al ). We therefore also 
compare our predictions with observations tha t estimate gas 
masse s based on the thermal emission of dust (jMagdis et al.l 
l2012al ) to circumvent this potential issue. Again, we find 
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Figure 3. Relation between halo mass and stellar mass of galaxies. (Left) Model predictions based on a linear £sfr — £h 2 relation, 
Equation J 1 3 1 ) - (Right) Model predictions based on the quadratic relation given by Equation 1141 . Solid lines are the model predictions 
for z = 0, z = 2, z = 5, and z = 8 (s ee legend). Symbols (c olor coded according to redshift, see legend) are empirical estimates based on 
the abundance matching technique (Bch roozi et alJ l2012aV Because the evolution of the Afhaio — M* relation is weak, we combine, for a 
given redshift z, the abundance matching results from the z — 1 (downward pointing triangles), z (squares), and z + 1 (upward pointing 
triangles) redshift bins. 



good agreement if the star formation - gas relation is linear, 
but stark disagreement if the relation is quadratic. 

The relation between stellar mass and halo mass of a 
galaxy is also sensitive to the functional form of the star 
formation - gas relation, see Fig. [3] Again, a linear star 
formation - gas relation leads to much better agreement with 
abundance matching results than a quadratic star formation 
- gas relation. The agreement is particularly good at high 
redshifts z > 4, where, as we show in the next section, the 
cosmic star formation history is sensitive to changes in the 
star formation - gas relation. 



3.2 The Cosmic Star Formation History 

In Fig. |4]we present the cosmic star formation history as pre- 
dicted by the model for a linear star formation - U~2 relation 
(Equation [T3J). Here, the black solid curve shows the contri- 
bution of all galaxies sufficiently bright (Mtjv,ab < —17.7) 
to be accessible to current observations at z ^ 4. We esti- 
mate the luminosity of a galaxy based on it s SFR and dust 
attenua tion. For the latter we use Table 6 of lBouwens et af] 
(|2012bh and we extrapolate the multiplicative dust extinc- 
tion to 1 (no dust extinction) at z = 10. For z < 3 we as- 
sume that the dust extinction remains constant at its value 
at z = 3. The red solid curve shows the co ntribution of all 
galax ies down to the photo-ionization limit (|Okamoto et al.l 
200gj) to the SFRD. The corresponding dashed curves show 
what happens if we assume that all gas in the disk is molecu- 
lar, i.e., for E = Eg in Equation (|13[) . T he figure also include s 
compilations o f observ ation al data by Oesch et al.l (|2012h . 
iBouwens et"afl l|2012bl ). and lBehroozi et all (|2012bT ). 



An important result of Fig. [3] is that a linear star for- 
mation - gas relation predicts a cosmic star formation his- 
tory that is in good agreement with observational data. In- 
terestingly, given the current observational limits the cos- 
mic star formation history at z > 2 changes only by a 
moderate amount depending on whether star formation 
is based on molecular hydrogen or not. We can rephrase 
this result, namely that quenching of star formation in low 
metallicity and hence H2-poor galaxies at high redshift is 
not very important for galaxies accessible to observations 
(Mtjv.ab < —17.7). It does, however, play a crucial role for 
fainter galaxies. In fact, considering all galaxies the cosmic 
star formation rate density at z ~ 10 is suppressed by 1.5 or- 
ders of magnitude if star formation occurs only in molecular 
gas compared with all the gas in the disk. Obviously, this will 
have implications for t he re-ionization by faint gal a xies a t 
such redshifts, see e.g., iKuhlen fc Faucher-Giguerel i|2012l ). 
which we plan to address in a follow-up study. 

In Fig. O we show the predicted UV luminosity func- 
tion and its evolution with redshift for a lin ear Esfr — Eh, 
relatio n and compare it with observations l|Bouwens et al.l 
l2012ah . The agreement is excellent considering the simplic- 
ity of the model and the various observational uncertain- 
ties. At z ~ 8 the observed luminosity function appears to 
have a somewhat larger normalization compared with our 
model predictions, while no such offset is seen at the other 
redshifts, which could indicate that the z ~ 8 observations 
are affected by cosmic variance. At low luminosities (not 
shown) the UV luminosity function turns over as a result of 
the suppression of star formation in low mass halos. Unfor- 
tunately, the magnitude at which the turn-over takes place 
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Figure 4. The cosmic star formation history as predicted by 
the model for a linear Esfr — Sh 2 relation (Equation 1131 with 
E = Sh 2 ; solid and dot-dashed lines) and a linear Ssfr — S g rela- 
tion (Equation 1131 with £ = S E ; dashed lines), r espectively, with 
a depletion time of 2.5 Gyr jBigiel et al.ll201lh . The H 2 abun- 
dance is computed either based on the current gas metallicity 
(solid lines) or by assuming that the gas has solar metallicity 
(dot-dashed lines). Black lines (lower three curves) show the con- 
tribution of galaxies bright enough to be accessible to current ob- 
servations, while red lines (upper three cu rves) include all galax - 
ies down to the photo-ionization limit ( |OkamotcT ct al. 2008). 
The contribution of galaxies below the current detection limit 
depends sensitively on whether star formation correlates prefer- 
entially with H2 or with the total gas mass in galaxies. Points 
with er ror bars are observational constraints from lBouwens et alj 
(2012b) and should be compared with the black lines. The dot- 
dashed line and the shaded region in the background show the 
mean and the typical error of the cosmic star form ation history, 
respectively, as estimated by Bchroozi et al. (2012b) using abun- 
dance matching and a compilation of observational data from the 
literature. Predictions based on a linear EgpR — £jj 2 relation with 
a ~ 2.5 Gyr H2 depletion time are in good agreement with the 
observed cosmic star formation history. 



depends sensitively on model parameters, e.g., on the disk 
sizes (parametrized by the spin parameter A), and is there- 
fore not a robust prediction of the model. For instance, at 
z ~ 4 the UV luminosity function is predicted to break at 
Mtjv ~ — 15.5 if we use a linear Esfr — £h 2 relation with 
our default choice of model parameters. If, however, we de- 
crease A from 0.07 to 0.05 the bright end of the luminosity 
function remains unaffected, while the turn-over shifts to 

Mtjv 12.5. 

In Fig. [6] we show plots analogous to Fig. [4] for four 
non-linear star formation - gas relation: a moderately sub- 
linear relation (Equati on 11511. the mod e rately super-linear 
relation postulated by iKrumholz et al.l l)2009bl) (Equation 
1 1611 . t he "universal relation" proposed by IKrumholz et 
((2012J) (Equation I17[) . and a quadratic relation (Equation 
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I14p . The figure clearly demonstrates that a super- linear star 



Figure 5. The UV luminosity function as predicted by the model 
for a linear EgpR — Ejj 2 relation with a 2.5 Gyr H2 depletion time. 
SFRs are converted into UV magnitudes according to Mtjv = 
-2.5 log 10 (SFR/ / dust [ik^yr" 1 ])- 18.16. A 0.1-125 M Q Salpeter 
initial mass function (IMF) and a constant SFR over the last 
> 100 Myr is assumed for the conversion. Solid lines correspond to 
UV magnitudes that include extinction by dust (/dust > l)i while 
dashed lines show the extinction free predictions (/d us t = !)■ We 
use the multiplicative du s t atten uation factor /d us t( 2 ) a s given in 
Fig. 14 of lBouwens et aL 1(20123) for z = 4 and z = 6 (the latter is 
assumed to hold for z ~ 6 — 9), and extend it to faint magnitudes 
such that /dust — 1- The hatched regions are measurements of 
the UV luminosity functions bv lBouwens et al] ll2012bl ). We use 
the Schechter function fits provided in their Table 3 and vary 
the fit parameters within their quoted error bars to create the 
hatched regions. The agreement between model predictions and 
observations is reasonable given the simplicity of the model, the 
uncertainties in the dust extinction corrections and / or changes in 
normalization due to cosmic variance. 



formation - gas relation results in SFRDs higher than ob- 
served at z ~ 4—10. The discrepancy is particularly large for 
strongly non-linear star formation-gas relations, such as the 
those given by Equations (|14|l and (|17[1 . These relations re- 
sult in SFRDs that are inconsistent with observations, e.g., 
at z = 8 (2 = 10) they are one (two) orders of magnitude 
above observations. Moderately non-linear relations such as 
those given by Equation (|16p and Equation II 7p are in ten- 
sion with observations, but are probably not excluded given 
the uncertainties of the model. 

Similar to Fig. |3J the observable SFRD shows no large 
changes depending on whether star formation takes place 
only in molecular gas or not. A sub-linear star formation 
- gas relation leads to the biggest difference (~ 0.5 dex), 
while there is almost no change if the relation is strongly 
super-linear. Furthermore, for a strongly super-linear (e.g., 
quadratic) star formation - gas relation there is little dif- 
ference in the star formation history even if galaxies below 
the observational limits are included. Hence, in this case, 
quenching of star formation in low metallicity and, thus, 
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Figure 6. Same as Fig. [3] but for non-linear SgFR — S g and Esfr — Sh 2 relations. (Top left) Predictio ns for a moderately sub - 
linear relation, see Equation (15) . (Top right) Predictions for the moderately super-linear relation suggested bvlKrumholz et alj (2009bh . 
see EquationJTH). (Bottom left) Predictions for the "universal" sta r formation - gas relatio n sug gested bv lKrumholz et alj d2012h . see 
Equation 11TI . (Bottom right) Predictions for a quadratic relation t Ostriker fc Shetty||201lh . see Equation H14II . The observable cosmic 
star formation history (black lines) is highly sensitive to the functional form of the star formation - gas relation. 



H2-poor galaxies is not affecting the cosmic star formation 
history. 

Before we address this last point, let us understand why 
the SFRD depends sensitively on the functional form of the 
star formation - gas relation if star formation is based on 
the abundance of H2. We show in Fig. the evolution of 
the typical gas surface density (here E g (r 1 / 2 ), the surface 
density at the gas half mass radius) of galaxies residing in 
halos of different mass. At z ~ 20 — 30 this surface density is 
typically in the range 300 — 1000 Mq pc -2 . As the universe 



expands, gas surface densities in galaxies decrease and reach 
values in agreement with observations of nearby galaxies (~ 
10 Mq pc ~ 2 for galaxies in ha los with masses ~ 10 11 — 10 12 
Mq, e.g.. lWong fc Blitz| [2002). Galaxies in halos above the 
virial shock quenching mass (> 2 x 10 12 Mq) experience a 
fast decline in their gas surface density at low redshift caused 
by the shut-down of gas accretion. 

We also include in the figure the evolution of the criti- 
cal surface density E g (/H 2 = 0.5) required to turn half the 
gas molecular. We note that there is a one-to-one correspon- 
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Figure 7. Evolution of the gas surface density at the disk half mass radius (E g (ri/ 2 ), solid lines) and evolution of the surface density at 
which the molecular fraction is 50% (E g (/n 2 = 0.5), dot-dashed lines) for halos of different masses (see legend). (Left) for a linear star 
formation — gas relation (Equation ll3l with E = Eh 2 ). (Right) for a quadratic star formation — gas relation (Equation 1141 with E = En 2 ). 
The expansion of the universe, the consumption of gas, and the quenching of gas accretion result in a reduction of E g (r 1/ / 2 ) with time. 
Star formation increases the metallicity and thus decreases E g (/n 2 = 0.5) (Equation [Sj. A quadratic star formation - gas relation leads 
to early enrichment (to ~ 0.1 Zq already at z > 10). Enrichment proceeds much more slowly in case of a linear Esfr — Eh 2 relation. 



dence between this critical surface density and the metallic- 
ity of the gas, see Equation ((H). Initially the critical sur- 
face density is large ~ 3000 Mq pc -2 (corresponding to 
■Zigm = 2 x 10~ 5 ), but it decreases with time as metal en- 
richment by star formation boosts the molecular content of 
the gas. Clearly, the choice of the star formation - gas re- 
lation has a large impact on the evolution of this critical 
surface density or, equivalently, the metallicity. If the star 
formation - gas relation is linear, halos with z = masses 
below ~ 3 x 10 10 M e are unable to significantly increase 
the metallicity of the gas. In contrast, efficient enrichment 
can take place in such halos at early times (z > 10) if the 
relation is non- linear. 

In Fig. [HJ we compare the ratio between the typical gas 
surface density, E g (r 1( / 2 ), and the critical surface density, 
^g(/n 2 = 0.5). H2-based star formation can only take place 
throughout the disk if this ratio is > 1. If the star forma- 
tion - gas relation is linear, gas in galaxies with z = 
halo masses below ~ 3 x 10 10 M Q , remains effectively un- 
enriched and thus atomic. Hence, star formation is almost 
completely quenched in halos of such mass. However, galax- 
ies in more massive halos are able to turn a substantial frac- 
tion of their gas reservoir into H 2 . In fact, the condition 
S g (ri/ 2 ) ~ E g (/H 2 = 0.5) is satisfied increasingly earlier for 
halos of higher mass, e.g., at z ~ 10 (z ~ 4.5) for halos 
with 2 = masses of 10 12 Mq (10 11 Mq). If the star for- 
mation - H 2 relation is quadratic, galaxies generally reach 
metallicities of Z ~ 2 x 10 -3 (a tenth solar) quickly after a 
single Pop III star initially enriches the gas to Z ~ 2 x 10 -5 . 
Hence, most galaxies, even those in low mass halos, have 
fully molecular gas disks and are able to form stars effi- 



ciently at early times. This also explains why the metallicity- 
dependent quenching of star formation has little impact on 
the cosmic star formation history in this case. Early enrich- 
ment (~ 0.1 Zq) and high gas surface densities (> 100 Mq 
pc -2 ) ensure that most of the gas in galaxies is molecular at 
z > 4, even in low mass halos, and thus it does not matter 
whether £h 2 or S g is the gas surface density relevant for 
star formation. 

Returning to the discussion of the linear star formation 
- gas relation, it is important to point out that the lack 
of early enrichment is not the reason that < 3 x 10 10 Mq 
halos are inefficient star formers at observationally accessible 
redshifts z < 10. To show this we re-run the model for the 
linear relation with a 100 times higher initial gas metallicity 
of the halos, i.e., with Z = 2 x 10" 3 instead of 2 x 10" 5 , 
see Fig. |8j Clearly, the ratio between typical and critical 
gas surface density is not affected in a significant way at 
z < 10. In particular, galaxies with z = halo masses below 
3 x 10 10 Mq still do not contribute to star formation at these 
redshifts. 

We conclude that while star formation follows the metal 
enrichment of the gas, the evolution of the latter is primarily 
driven by the functional form of the star formation - gas 
relation. A non-linear relation leads to fast enrichment and, 
hence, to galaxy disks that are fully molecular and have 
large star formation rates even in low mass halos at high 
redshifts. In contrast, a linear relation delays the enrichment 
of the gas and therefore limits star formation to sufficiently 
massive halos. 

This leaves us with one open question, namely how can 
the star formation - gas relation have such a strong impact 
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Figure 8. Evolution of the ratio between the typical gas surface density in a disk, Sg(ri/2)i an d the surface density at which the 
molecular fraction is 50%, £ g (/H 2 = 0.5), see Fig. [7] A ratio significantly below (above) unity means that most of the gas in the disk of a 
galaxy is atomic (molecular). (Left) for a linear star formation - gas relation (Equation 1131 with £ = £h 2 )- G as surface densities become 
sufficiently large to form molecules and stars only in massive halos at sufficiently late times. (Middle) for a quadratic star formation - 
gas relation (Equation 1141 with E = £h 2 )- A strongly super-linear star formation — gas relation allows the formation of fully molecular 
gas disks at earlier times and in halos of lower masses and, hence, boosts the star formation rate at the relevant, and observationally 
accessible, redshift range z ~ 4 — 10. (Right) same as left panel, but with the metallicity of the gas in the galaxy disks initially set to 
Z = 0.1 Zq. The 2 < 10 evolution is hardly affected by this dramatic change, demonstrating that the lack of early enrichment is not the 
reason that a linear star formation - gas relation is able to reproduce the observed cosmic star formation history at z ^ 10. 



on the metal enrichment in halos? The metallicity evolu- 
tion is determined by the SFR (increases metal mass and 
decreases gas mass; Equation 1 12|) and the gas accretion rate 
(increases gas mass and metal mass; Equation[4]|. Rearrang- 
ing these equations it is easy to show (see appendix [X} that 
the metallicity Z is constant, i.e., Z = 0, if 

_ M*,fo Tm _ Z — ZlGM 

r = ^w~7~y(i-^)(i-<r ( ' 

Given a SFR to gas accretion rate ratio r, the equilibrium 
metallicity is therefore 

Z cq {r) = ry{l-R)(l-() + Z IG M. (21) 

We demonstrate in appendix[X]that this equilibrium so- 
lution is (linearly) stable under most circumstances. Hence, 
for a given ratio r the metallicity of the ISM will increase or 
decrease until it reaches Z eq (r). In other words, the compe- 
tition between SFR and gas accretion rate drives the metal- 
licity evolution in galaxies. 

Equation (|21|l has a few interesting properties. First, it 
demonstrates that galactic outflows of warm/cold ISM ma- 
terial (parametrized by e ou t) play no role for the ISM metal- 
licity. The reason is that while outflows carry metals away 
they also remove a corresponding amount of gas leaving the 
equilibrium metallicity unchanged. But since the metallicity 
of the ISM is likely close to Zeq(r), it also remains unaf- 
fected. Second, Z eq (r) depends only weakly on halo mass 
via £ and, assuming Zigm = const, not on redshift. Hence, 
an analysis of the mass or redshift dependence of the metal- 
licities of high redshift galaxies (combined with measured 
SFRs and mass accretion rates, estimated, e.g., from their 
halo masses) could be used to learn more about supernova 
driven outflows (parametrized by £) or the IGM enrichment 
(•Zigm). 

In the next section, we study the connection between 



the star formation - gas relation and the driver of the metal- 
licity evolution in galaxies, the ratio of SFR to gas accretion 
rate. 



3.3 Gas Accretion and Star formation 

Previously we have demonstrated that the functional form 
of the star formation - gas relation has a profound impact on 
the SFR of galaxies at high redshifts and we have attributed 
this effect to the altered metal enrichment. In this section, 
we will show that the metal enrichment is intrinsically linked 
to the star formation - gas relation via the balance between 
star formation and gas accretion in halos. 

In Fig. [5] we plot the evolution of the ratio of the accre- 
tion time (gas mass divided by the gas accretion rate) to the 
gas depletion time (gas mass divided by SFR), i.e., the ratio 
of the SFR and the gas accretion rate, for halos of a range 
of masses. At z > 3, a linear star formation - gas relation 
(Equationsll3p results in a gas depletion time that is signif- 
icantly longer (it is 2.5 Gyr if all gas is molecular and even 
longer otherwise) than the gas accretion time (~ 2-3x 10 7 
yr at z = 30, ~ 1 - 2 x 10 s yr at z = 10, 2 - 7 x 10 s 
yr at z = 5). Hence, star formation is "depletion limited" 
at those redshifts. Changing the functional form of the star 
formation gas - relation alters the gas depletion times, but 
has only a small effect on the gas accretion rates, which are 
determined primarily by the mass accretion rates onto dark 
matter halos. A quadratic star formation gas - relation, for 
instance, results in gas depletion times as low as ~ 2 x 10 8 
yr at z > 10. As a consequence, star formation is "accretion 
limited" (the SFR is of the same order as the gas accretion 
rate) over most of cosmic history in this case. 

Equation (120[) allows us to understand the consequences 
of a depletion limited vs an accretion limited star formation. 
Consider a gas disk with a given metallicity Z. Such a gas 
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Figure 9. Ratio between gas accretion time (gas mass divided by gas accretion rate) and gas depletion time (gas mass divided by SFR), 
or equivalently the ratio between SFR and gas accretion rate, for galaxies in halos with present day masses in the range 10 10 - 10 13 M 
(see legend). Egppt is assumed to scale either linearly (left panel) or quadratically (right panel) with Sh 2 - Star formation is limited by 
gas accretion if the ratio is of the order of unity. The ratio gets very large in high mass halos at low redshift due to virial shock quenching 
of gas accretion, see Section [21 ■ The dot-dashed lines shows the ratio of accretion to depletion time that is required to keep the gas 
metallicity constant at the value indicated in the figure. These lines are computed using Equation (1 20 P with £ = 0.9 for Z = 0.1 — 0.01 
Zq and C = 0.5 for Z = Zq. A ratio of accretion to depletion time above a given dot-dashed line with label Z indicates that metal 
enrichment via star formation outmatches the accretion of gas of low metallicity and, as a net effect, that the metallicity can raise above 
Z. A linear star formation - gas relation (Equation ll3|l has a gas depletion time that is significantly larger than the gas accretion time at 
z > 3. Metal enrichment and, hence, star formation is delayed until such times. In contrast, a quadratic star formation - gas relation has 
a depletion time that is short enough for star formation to catch up with the accretion rate early on. This leads to a fast and sustained 
enrichment to ~ 0.1 Zq at high redshift and, consequently, to high molecular fractions and substantial star formation over much of 
cosmic history. 



disk will increase or decrease its metallicity until it matches 
the metallicity Z eq specified by Equation (|2ip for the given 
SFR to gas accretion rate ratio r. For a stellar mass loss 
fraction of R = 0.46, a stellar metal yield of y = 0.069, 
and £ = 0.9 (true for all but the most massive halos at 
z ~ 4 — 10), we obtain Z cq = 0.1 Zq for r — 0.53, but 
Z — 0.01 Zq for r = 0.048. Therefore, if star formation is 
accretion limited the gas will get enriched to ~ 0.1 Zq or 
more on the gas accretion time scale. A rapid enrichment to 
Z ~ 0.1 Zq and accretion limited star formation has indeed 
been seen in many hydrodynamical simu lations with a super - 
linear star formation model, see, e.g., iDave et al.l (|2006ri ; 
iFinlator et al.l (|201ll ). The equilibrium metallicity is small 
if r is small, i.e., if the depletion time is much longer than 
the gas accretion time. Hence, metal enrichment is delayed 
until z ~ 2 — 3 in case of a linear star formation - gas relation 
with a sufficiently large (>Gyr) depletion time. 

To summarize, the functional form of the star formation 
- gas relation determines whether star formation in a given 
halo at a given time is depletion limited or accretion limited. 
This is crucial, because the evolution of the metallicity in a 
halo is driven by the competition between gas accretion and 
gas depletion (Equation I21|) . A linear star formation - gas 
relation delays metal enrichment until late times and, hence, 
suppresses star formation at high redshifts. In contrast, a 



strongly super-linear star formation - gas relation boosts 
metal enrichment and star formation at high redshifts. 

A quantity similar to the ratio between SFR and gas ac- 
cretion rate has been studied empirically by Be hroozi et al.l 
(|2012al ) using abundance matching techniques. In their work 
the baryon equivalent accretion rate (BEAR; /bMhalo) is 
used as a proxy fo r the unknown gas ac cretion rate. They 
find (see Fig. 1 in iBehroozi et al J 120123 ) that the SFR to 
BEAR ratio is roughly constant out to z ~ 3 and specu- 
late that it might be a time-independent for much of cosmic 
history. 

In Fig. [10] we show our predictions for the SFR to BEAR 
ra tio using a similar color scheme and axis ranges as in Fig. 1 
of IBehroozi et al.l (|2012al ). We compute baryon equivalent 
accretion rates based on t he halo mass evolution formula 
provided in the appendix of lBehroozi et all (|2012al ). We find 
that a linear, tb-based star formation gas - relation can 
reproduce their Fig. 1 reasonably well. Most importantly, 
the ratio between SFR and BEAR peaks around halo masses 
of ~ 10 Mq at z < 2. Some features seen in Fig. 1101 such as 
the complete lack of star formation in sufficiently low mass 
halos, may be related to the fact that the model follows the 
average evolution of galaxy properties and, at least in its 
current form, ignores the possibility of scatter caused by, 
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Figure 10. Ratio between SFR and baryon equivalent accretion 
rate (/b^Whaioj BEAR) of galaxies as function of cosmic time and 
halo mass assuming a linear Ssfr — ^H 2 relation. For z < 2 
the ratio peaks at a halo mass of appr oximately 10 11 Mp, i n 
reasonable agreement with observations l|Behroozi et alj|2012ah . 
Galaxies in low mass halos (A4ai < 2 X 10 10 M Q at z ~ 0) 
have their star formation quenched due the scarcity of molecular 
hydrogen in such halos. In addition, halos stop accreting gas once 
they become massive enough (the upper dashed line) to support 
a virial shock that preve nts cold gas accretion onto the galaxy 
l|Dekel fc Birnboim|[200r3) . At z < I star formation is confined to 
a narrow range of halo masses of approximately a few times 10 10 
M© to 10 12 M Q l|Bouche et alj|20ld: [Krumholz fc Dekelll2012h . 



e.g., the variability of the accretion rates or scatter in the 
disk sizes of galaxies. 



3.4 An equilibrium view on galaxy evolution 

Not only does the ratio between SFR and gas accretion rate 
regulate the metallicity in galaxies (see Equation [21]), it also 
determines the gas fraction f s — M g /M^ a i a , stellar fraction 
/ s = M*/Mh a io, and the metal fraction fz = Mz/Mhaio of 
galaxies. 

Specifically, we can show using Equations [2] [3] [4] and [5] 



that for a given ratio r 



M t,brm 



fg = if f s = /g,eq(r) with 

/g,eq(r) = [1 - r(l - R + e out )] em/g,in/b, 

f s = if f s = /s, cq (r) with 

/s,eq(f) = (1 - ,fg,in)/b + £in/g,in/b(l - R)r, 

fz = if fz = fz,cq(r) with 

/z.oqO) = Zeq(r)/g,eq(r), . 



(22) 
(23) 
(24) 



The equilibrium values of other properties, such as the gas- 
to-stellar mass ratio / gs = f s /f s or the baron fraction f s +f a , 
can be obtained by combining these equations. For instance, 



/gs,eq( r ) 



,fg,cq(r) _ l-(l-R + e 

/s,eq(r) 



+ (1 - R)r 



(25) 



We show in appendix|X]that the equilibrium state spec- 
ified by Equations H22l) - (|24| ) is (linearly) stable under most 
circumstances for SFR to gas accretion rate ratios less than 
1/(1 — R + e ut)- If, however, the galaxy is in a condition in 
which the H2 fraction is extremely sensitive to small changes 
of the metal abundance of the ISM and star formation de- 
pends on the presence of molecular hydrogen, a run-away 
increase or decrease in H2 abundance might take place. The 
latter ceases, and the galaxy reaches a new equilibrium, as 
soon as the H2 fraction becomes sufficiently insensitive to 
changes in metal abundance. There can be no equilibrium 
state for r > 1/(1 — 7? + e out ), because it would lead to neg- 
ative equilibrium gas fractions, see Equation 1221 However, 
r > 1/(1 — R+ e ou t ) is a physical possibility and applies, e.g. , 
to massive galaxies at low redshifts that have their gas ac- 
cretion quenched, but are still star forming. These galaxies 
are obviously not in an equilibrium state, but rely on their 
diminishing gas reservoir for star formation. 

The stability of the equilibrium state implies, for in- 
stance, that if / g < f s ,eq(r) (and it is sufficiently close to the 
equilibrium value), the gas fraction increases and continues 
to do so until it reaches f g ,eq(r). Similarly, if / g > / g ,eq(r), 
the gas fraction decreases until it matches the equilibrium 
value. We can therefore draw important conclusions for the 
gas fractions and gas-to-stellar mass ratios from Fig. [5] A 
small SFR to gas accretion rate ratio of, e.g., r = 0.048 cor- 
responds to an equilibrium metallicity of Z — 0.01 Zq, to 
an equilibrium gas fraction / g ,eq//b ~ 0.83 (0.46) and to 
an equilibrium gas-to-stellar mass ratio / ga ,eq ~ 6.8 (0.9). 
Here we assumqj / g ,m = 0.9 (/ g ,m = 0.5) and the default 
parameters of the model: ei n = 1, R = 0.46, e ou t = 1. Hence, 
depletion limited star formation implies low metallicities, 
large gas fractions, and large gas-to-stellar mass ratios. In 
contrast, equilibrium metallicities are higher and equilib- 
rium gas fractions and gas-to-stellar mass ratios are much 
lower if star formation is accretion limited. For instances, 
a SFR to gas accretion rate ratio r = 0.53 corresponds to 
Z cq = 0.1 Z Q , /g, cq //b ~ 0.17 (0.09) and / gs , cq ~ 0.46 (0.14) 
for/ g> = 0.9 (/ g , in = 0.5). 

In the context of our model the baryonic properties of 
galaxies (of given halo mass and at a given redshift) are 
specified by fz, f a , and / s . In more abstract terms, the set 
(fz, f 3 , fs) represents the baryonic state of a galaxy. Equa- 
tions (|22[) - (|24[) denote a curvqj of fixed points parametrized 
by r in this baryonic state space. We demonstrate in the ap- 
pendix that if {fz, f g , fs) is reasonably close to one of the 

4 The parameter / g ,i n can affect the predicted equilibrium gas 
fractions and gas-to-stellar mass ratios noticeably. The values 
chosen here are for illustration purposes, but we stress that we 
compute / gl i n self-consistently in the model (see Section[2]l. 

5 This is only approximately true, because some of the terms in 
the Equations J22D - (1241 i dependent on more than just r, e.g., 
on the halo mass. This also implies that different galaxies will in 
general follow somewhat different curves of fixed points during 
their evolution, see Fig. 1111 
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fixed points it will approach it. As the cosmic accretion rate 
onto a halo changes, r changes as well and the galaxy transits 
smoothly from one equilibrium state to the next. Strong per- 
turbations, such as major mergers, will also allow a galaxy 
to leave its current fixed point, but we would expect that 
such events are relatively rare. Hence, galaxy evolution is 
largely an evolution along this curve of fixed points driven 
by the external modulation of the gas accretion rate. 

This clearly idealized picture assumes that galaxies can 
adapt on short timescales to changes in the SFR to gas ac- 
cretion rate ratio r. The matter accretion time iWhalo/-Whalo 
is a natural timescale for gas fraction, stellar fractions, and 
metallicities to reach their equilibrium values, but the pre- 
cise value depends on the parameters of the model, e.g., on 
the exponent of the star formation - gas relation, see ap- 
pendix [X] Given that the accretion time is of the order of 
the Hubble time at a given redshift, we expect that the stel- 
lar fractions, gas fractions and metallicities lag somewhat 
behind their equilibrium predictions, i.e., that gas fractions 
are somewhat underestimated by the equilibrium formalism, 
while stellar fractions and metallicities are overestimated. 
Fig. [11] confirms this point. Here, we show how halos evolve 
at redshift z > 2 in the f g -Z projection of the baryonic 
state space. We plot both the metallicities and gas frac- 
tion as computed by the model, and also the corresponding 
equilibrium quantities based on Equations (|2ip and (|22p . 
The equilibrium formalism describes the overall trend be- 
tween / g — Z well, but for a given r it underpredicts the 
gas fractions and overpredicts the metallicities significantly. 
While properly correcting Equations (12211 - (125[) for the time 
lag is beyond the scope of this paper, we point out that 
the disagreement can be reduced if r is scaled by a factor 
~ 0.45[2- (2-n) 2 ] before applying Equations ([21} and ([22)) . 
see Fig. 1111 Here, n £ [1,2] is the power- law exponent n of 
the star formation - gas relation. 

We note that our equilibrium picture of galaxy forma- 
tion differs fr om the "equilibrium condition" of a s teady 
state model (iFinlator fc Pavel |200S| ; iBouche et ail l20ld ; 
button et aLll2010l ; lDave et al-lboil Tin which the SFR in a 
given halo is balanced by the net gas accretion rate. Most 
importantly, we demonstrate that galaxies can be in a sta- 
ble equilibrium even if the "equilibrium condition" is not 
even remotely satisfied. Specifically, galaxies can have little 
star formation and can massively build up their gas reser- 
voirs and, yet, be in a proper equilibrium state. The rea- 
son is that the "equilibrium condition" of the steady state 
model is equivalent to the assumption that the gas mass in 
a galaxy remains constant. In contrast, we show that the 
actually stable equilibrium condition is that the gas frac- 
tion remains constant. The special equilibrium state with 
r = 1/(1 — R + e ut) satisfies the "equilibrium condition" of 
the steady state model. Indeed a similar relation for r is of- 
ten assumed to hold generally in t he ste ady state approach, 
see e.g., Equation 2 of Dave et alj (|2012l ). However, as Equa- 
tion (|22|l shows the equilibrium gas fraction, and, hence, gas 
mass, in this equilibrium state is zero, and, hence, gas-rich, 
star forming galaxies at high redshift will not be in this par- 
ticular state. 

The equilibrium view presented in this section shows 
why metallicities, gas fractions, and stellar fractions of 
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Figure 11. High redshift evolution of gas fractions and metal- 
licites of halos with z = masses of 10 11 Mq (red lines), 10 12 
Mq (black lines), and 10 13 Mq (blue lines), respectively. The 
solid lines and stars show the evolution as predicted by the 
model presented in Section [2] for a linear Ssfr — ^H 2 relation. 
The dot-dashed lines and the crosses (dashed lines and squares) 
show the equilibrium predictions of Equations II21I I and i ]22[ ) with 
r = M» >form /Af gi i n (r = 0.45M» iform /M gi i n ). Symbols show pre- 
dictions at redshifts 10, 8, 6, 4, 3, and 2 (from top left to bottom 
right). The equilibrium ansatz predicts the right locus in the / g 
vs Z plane, but for a given r underestimates / g and overestimates 
Z. This disagreement is an indicator that galaxies have not suf- 
ficient time to fully reach the equilibrium solution. The simple 
rescaling of r by a factor ~ 0.4 — 0.5 improves the accuracy of the 
equilibrium formalism at z ^ 2. 



galaxies are strongly correlated. They are all coupled via 
a single parameter - the ratio between the SFR and the gas 
accretion rate, see Equations (|22[) - (|24[) . Since the SFR de- 
pends on these galaxy properties, it follows that for a fixed 
functional form of the star formation - gas relation the evo- 
lution of global galaxy properties is regulated by just the 
gas accretion rate. Galactic outflows modify the equilibrium 
gas fractions (see Equation 122 fl , but the primary regulator 
of star formation and other global galaxy properties is the 
gas accretion rate. 



3.5 Caveats 

Our findings come with a number of potential caveats. 
First, our approach is based on the assumption that a sin- 
gle, unchanging star formation - gas relation, potentially 
parametrized by local galaxy properties such as the gas 
metallicity, or the free-fall time, etc., holds from early times 
until the present. Second, the model obviously idealizes 
many aspects of star formation in high redshift galaxies. 
For instance, it assumes that the gas is always arranged in 
an exponential disk (likely not true during major mergers) 
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Figure 12. Predictions of an alternative model in which the mass ejection parameter e ou t increases at z > 2. The star formation - gas 
relation is assumed to be quadratic (Equation [14} . The panels show from left to right the cosmic star formation history, the gas-to-stellar 
mass ratio vs stellar mass, and the relation between stellar mass and halo mass. The lines and symbols are as in Fig. [4] Fig. [2] and Fig. [3] 
for the left, middle, and right panel, respectively. The ad-hoc scaling of e ut with redshift cures the excess of SFRDs at high z otherwise 
predicted for a quadratic star formation - gas relation, see Fig. [6] However, it does not help in bringing the M^aio — relation or the 
gas-to-stellar mass fractions in agreement with observations. 



and that accreted halos donate all their cold gas to the cen- 
tral galaxy. Furthermore, since the model only makes pre- 
dictions for the mean evolution of halos of a given mass, 
it ignores variations in the accretion histories and possi- 
ble environmental dependences. Such, admittedly interest- 
ing, complexities are best studied in numerical simulations. 
However, demanding requi rements in terms o f simulation 
volume and resolution (see lSchave et ai1l2010l for a resolu- 
tion study) limit the systematic exploration of the parame- 
ter space of feedback and star formation models using this 
approach. Nonetheless, we plan to compare the predictions 
of our model with dedicated high resolution cosmological 
simulations in future work. 
Additional caveats are: 

• Model parameters: We u se the same set of parameters as 
in lKrumholz fc D ckel (2012). Most of these parameters have 
a theoretical or empirica l justification (we refer to reader to 
iKrumholz fc Dekelll2012T ). When we vary the model param- 
eters within their justifiable ranges, we find that the mass 
ejection rate into the IGM per SFR, e out , the metallicity of 
the IGM, Zigm, and the spin of the halo, A, have the largest 
effect on the predicted cosmic star formation history. The 
impact of these changes is largest if a linear star formation - 
gas relation is used, and smaller if the relation is highly non- 
linear and, hence, star formation is accretion limited early 
on. However, even for a linear relation we can vary each of 
these parameters by a factor of ~ 3 without any qualitative 
changes in the model predictions. For instance, increasing 
tout by such a factor decreases the SFRD at z ~ 1 — 2 by 
~ 0.5 dex, but leads to little change at higher redshifts (e.g., 
at z — 4 the decrease in the SFRD is only 0.2 dex). Reduc- 
ing Zigm by a factor 3 delays star formation in low mass 
halos and results in SFRDs that are reduced by a similar 
factor. Decreasing A from 0.07 to 0.05 increases the SFRD 
at z ~ 4 — 10 by a factor ~ 2 — 3. The predicted observ- 
able cosmic star formation history is less affected (< 0.2 dex 
typically). 



• Missing physics: Given the simplicity of the model it 
is possible that we are missing an important physical pro- 
cess that would affect the model predictions and potentially 
change our conclusions. For instance, a mechanism that in- 
creases the mass ejection rate of the winds with redshift 
would help to mitigate some of the problems of a highly non- 
linear star formation - gas relation. We demonstrate this in 
Fig. [I2]where we show the model predictions for a hypothet- 
ical scenario in which e out = 1 for z < 2 and e out = 2 — 1 for 
z ^ 2. Clearly, by choosing the scaling of e out properly even a 
quadratic star formation - gas relation can be brought into 
agreement with observations of the cosmic star formation 
history. However, this does not necessarily mean that other 
observational constraints are met as well. The figure shows, 
for instance, that this particular model fails in reproducing 
the Mhaio — M* relation at any redshift and the gas-to-stellar 
mass ratios at z > 1. In contrast, a simple, linear star for- 
mation - gas relation with an H2 depletion time of 2.5 Gyr 
(as observed in the local universe) does not require the as- 
sumption of additional ad-hoc (and fine-tuned) physics to 
reproduce global galaxy properties from high redshifts until 
today, see Section [3. II 

• Metal mixing: We assume that metals in the disks of 
galaxies are efficiently mixed, i.e., that the gas has a spatially 
uniform metallicity. In reality, we expect some metallicity 
gradients, since accretion of metal-poor gas occurs prefer- 
en tially at the outer ed ge of the disk. Numerical simulation 
bv lKuhlen et al l (|2012h indicate that the metallicity differ- 
ence is only modest for the majority of the gas in the disks 
of high redshift galaxies. In order to test the importance 
of metal mixing we modified the model to allow only local 
enrichment (i.e., metals produced in a given radius bin are 
added only to that radial bin) . This somewhat extreme case 
of no-mixing lowers the predicted SFRDs by ~ 3 if we use 
a linear star formation - gas relation, but makes almost no 
difference if the relation is strongly non-linear. 

• Dust: The model assumes that the dust-to-gas ratio 
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(which determines the H2 abundance) scales proportional 
to the metallicity. At low metallicities (Z < O.IZq) it 
is unclear whether this assumption is justified. A poten- 
tial concern is that at high z the dust production may 
lag behind the metal production, limiting the formation of 
H2, and thus star formation. Empirically, there is evidence 
for a substantial presence of dust in high redshift galax- 
ies from absorption line studies of damped Ly-q systems at 
z ~ 1.5 — 3.5 (e.g., iPettini et al. I ll994l ; iProchaska fc Wolfe! 
2002) and from t he thermal emission of dust in quasar host s 
at z ~ 4-6 (e.g.. lOmont et al.ll200ll;jBertoldi fc Coxll2002h . 
Both AGB stars (e.g.. [Valiante et al-booi T and supernovae 
(e.g. , iTodini fc FerraralbOOll ) have been proposed as forma- 
tion channels for dust at such high redshifts. If supernovae 
dominate the dust production, we expect a dust production 
time scale of ~ 10 7 yr, which is shorter than the gas accre- 
tion time (~ 1 — 2 x 10 s at z — 10) and, hence, not a limiting 
factor. 

• Molecular hydrogen: The H2 abundances are e s timate d 
based on the equilibrium model bv lKrumholz et all (|2009al ). 
For Milky- Way like, turbulent ISM cond itions the H2 forma- 
tion times tH 2 ,form are short (~ 10 6 vr. Idover fc Mac Lowl 
120071 ). but they could be higher in high redshift galaxies. 
With tn 2 .form OC (Zp)- 1 ijHollenbach et al.lll97lf ) and gas 
densities in the high z galaxies that are likely an order of 
magnitude higher than in the Milky Way, we expect forma- 
tion times of < 10 8 yr for gas with metallicities Z ^ 2 x 10 -5 , 
which is sufficiently short to not limit star formation. In ad- 
dition, numerical simulations of star formation in the ISM 
show that decreasing the metallicity from solar to a hun- 
dredth sola r postpones the onset of star formation only by 
a few Myr jGlover fc Clarldl2012h . Hence, we do not expect 
that star formation is limited by the H2 formation time even 
in case of a highly non-linear star formation - gas relation. 



4 SUMMARY AND DISCUSSION 

We have studied the implications of the functional form of 
the star formation - gas relation for the cosmic star forma- 
tion history, for gas accretion in galaxies, and for the evolu- 
tion of global galaxy properties using a model that follows 
the average evolution of galaxies in growing dark matter ha- 
los across cosmic time (|Krumholz fc De kcl 2012). We specif- 
ically contra sted the prediction s for a linear star formation - 
gas rel ation dBigiel et al.[|2 011) with those for a quadratic re- 
lation |Ostriker fc Shettvll201ll 'l. but also t ested a few other 
functional fo r ms suggested in the literature (|Krumholz et al.l 
l2009bl . 120121 ; IShettv et aI1l2012l 'l. Our main findings are as 
follows: 

• We show that a linear Esfr — Eh 2 relation with a ~ 2.5 
Gyr depletion time results in a mass - metallicity relation, 
in gas-to-stellar mass ratios, and in a stellar-to-halo mass 
relation that is in reasonable agreement with observations. 

• The functional form of the star formation - gas relation 
strongly affects the cosmic star formation history at z > 4. A 
linear star formation - H2 relation suppresses star formation 
in halos with z — masses below ~ 10 11 Mq. If, however, 
the star formation - gas relation is highly non-linear, star 
formation can occur in halos with z = masses even below 



~ 10 10 Mq. As a result, SFRDs can differ by factors of ~ 30 
at z — 10 (~ 5 at z — 6), see Fig.|4]and[6l between the linear 
and the quadratic star formation - H2 relation. 

• The star formation - gas relation has a similar effect on 
the observable cosmic star formation history, i.e., the SFRD 
of galaxies with magnitudes acces sible to observations (h ere 
assumed to be Muv.ab < —17.7, iBouwens et afll2012b[ ). A 
linear star - H2 formation relation results in an observable 
cosmic star formation history compatible with observational 
data, while a quadratic relation strongly overpredicts the 
observable SFRDs at high redshift. 

• Whether star formation is based on H2 or on total gas 
leaves little imprint on the observable cosmic star formation 
history. There is also little difference in the full (including all 
galaxies) cosmic star formation history if the star formation 

- gas relation is highly super-linear. If the relation is linear, 
however, the choice of Hb-based vs total gas based star for- 
mation affects the SFRD at high redshift, e.g, by more than 
an order of magnitude at z = 8. 

• For the adopted linear Esfr — Eh 2 relation (Equation 
I13[) . the gas depletion time exceeds the gas accretion time 
(~ free-fall time of the halo) at high 2. During this time 
the SFRs are orders of magnitude smaller than the gas ac- 
cretion rates. Star formation becomes accretion limited, i.e., 
SFRs ~ gas accretion rates, only at relatively recent epochs 
(z < 2 — 3). We show that these predictions are in rea- 
sonable agreement with abundance matching predictions by 
iBehroozi et~al1 1201231 ). In contrast, star formation is close 
to accretion limited already at very high redshifts if the star 
formation - gas relation is quadratic. 

• The evolution of metallicities, gas fractions, and stellar 
fractions are driven by the ratio r between the SFR and gas 
accretion rate. In fact, given a star formation - gas rela- 
tion galaxies evolve towards a (linearly stable) equilibrium 
state with gas fraction, stellar fraction and metallicity de- 
termined only by the gas accretion rate. The equilibrium 
concept presented in this work differs from the steady state 
model suggested in the literature fe.g.. |Pave et alj2012t ) and 
shows that galaxies can be in equilibrium even if their gas 
masses are growing quickly. We argue that galaxy evolution 
can be interpreted as an evolution along a set of equilibria 
driven primarily by the change in the cosmic accretion rate. 

As pointed out above we find that the observed cosmic 
star formation history is well reproduced if the star forma- 
tion - gas relation is linear and has an H2 depletion time 
of ~ 2.5 Gyr. Given the simplicity of the model and the 
uncertainties in the observational data we do not rule out 
the possibility that the relation is actually slightly super- or 
sub-linear (e.g., has an exponent in the range ~ 0.85 — 1.3) 
or has a somewhat different depletion time. However, our 
findings are difficult to reconcile with the idea of a highly 
super-linear, especially quadratic, star formation - gas rela- 
tion. 

Many of our findings for a linear star formation 

- gas relation are almo st identical to those found by 
iKrumholz fc Dekell l|2012l ) for a non-linear star formation 

- gas relation, similar to Equation 1161 with an exponent of 
1.33 at high gas surface densities. However, their star forma- 
tion - gas relation uses a very large transition surface density 
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(E g = 850 Mq pc~ 2 compared with E g = 85 Mq pc -2 of 
Equation I16[) at which the scaling changes from sub-linear 
to super-linear. This matters because gas surface densities 
of high redshift galaxies are of the order of a few 100 Mq 
pc~ 2 , see Fig. [7] Hence, a transition surface density as large 
as Eg = 850 Mq pc -2 implies that a large fraction of the gas 
in galaxy disks is subject to a sub- linear star formation - gas 
relation and that Esfr scales super-linearly with Eh 2 only 
in the central regions (where surface densities are larger by a 
factor e 1,7 ~ 5.5 than the surface densities at the half-mass 
radius). Averaging over sub- and super- linear regions results 
in a close-to-linear relation between H2 masses and SFRs in 
galaxies, hence, to similar predictions as for the strictly lin- 
ear star formation - gas relation given by Equation (|13[) . 

iKrumholz fc Dekell |2012l ) argue that star formation 
based on H2 is an essential ingredient to bring the predicted 
SFRDs at high redshifts in agreement with observations. We 
find that metal-dependent quenching of star formation based 
on H2 has only a small impact (0.3 dex at z ~ 9 and less 
at lower redshift) on the observable SFRDs at high redshift, 
because it affects primarily faint galaxies below the current 
detection limits. However, metal-dependent quenching plays 
an important role for the total (i.e., including galaxies below 
the detection limits) SFRD if the star formation - gas re- 
lation is not highly super -linear. This is in agreement with 
IKrumholz k, Dekell l|2012ft who effectively use an approxi- 
mately linear star formation - gas relation (see above) to 
study the total SFRD. In case of a strongly super-linear 
relation, however, star formation in galaxies is almost al- 
ways close to being accretion limited even at 2 ~ 10. There- 
fore, galaxies get enriched early on and SFRs proceed at a 
rate that is effectively identical (to within a factor of 2 or 
so) to the gas accretion rates, independently on whether a 
Esfr — Eh 2 or a Esfr — E g relation is used. 

ISchave et al.l j2010) analyze the star formation history 
of the universe with large scale cosmological simulations and 
find little evidence for a change of the (total) cosmic star for- 
mation rate with metallicity. The star formation in their sim- 
ulations is modeled via a moderately non-linear "Schmidt 
law", i.e., a star formation - gas relation based on volumet- 
ric densities. They allow star formation only in regions with 
densities above a threshold that is either fixed (the default 
model) or scales with metallicity. However, their threshold 
is much lower than the d ensity required for H2 formation 
dGnedin fc Kravtsovllioill) a nd the metallicity scaling much 
weaker (oc log Z~\ ISchave! |2004 ) than the one expected 
to regulate the presence of H2 (approximately oc Z~ x , see 
Equations [8] [TT}. It is therefore difficult to properly judge 
the implications of their findings in the context of H2 based 
star formation. 

iKuhlen et al.l ([2012) study the evolution of the cosmic 
star formation history in simulations that follow the abun- 
dance of H2 directly (using Equations f8l 41 1 p and compare 
these to simulations that use a fixed density threshold for 
star formation. The chosen threshold (50 cm -3 ) corresponds 
approximately to the transition from atomic to molecular 
phase for gas of solar metallicity. The Schmidt law used in 
the simulations is moderately non-linear. Hence, we would 
expect that the two sets of simulations result in very similar 
predictions for the SFRs of galaxies bright enough to be ob- 



servable. Indeed, IKuhlen et all (|2012T l find that the SFRDs 
of galaxies with stellar masses larger than 10 7 ' 75 Mq change 
by (only) a factor of ~ 2. The differences are even smaller 
when galaxies are selected based on their UV luminosity 
(Muv.ab < —17.7; M. Kuhlen private communication). 

Another interesting question is the importance of stel- 
lar feedback f or the regulation of star formation at high z. 
ISchave et ail (|2010h found that increasing the exponent of 
the star formation - gas relation or decreasing the depletion 
time has little impact on the cosmic star formation history 
at z < 6 and the authors use this fact to infer that galaxies 
are able to regulate their star formation via feedback pro- 
cesses. However, their default star formation - gas relation 
is non-linear and, hence, increasing its exponent or decreas- 
ing its depletion time further has little impact if SFRs are 
already close to the accretion limit. We also argue, based 
on the equilibrium view of galaxy evolution presented in 
Section 13.41 against the interpretation that star formation 
is regulated by large scale galactic outflows, because for a 
given star formation - gas relation the gas accretion rate (or 
the matter accretion rate onto the halo) is the main regula- 
tory mechanism of star formation. Of course, star formation 
may be feedback-regulated within the ISM. In particular, 
the existence of a star formation - gas relation that holds 
on ~ kpc scales may be a consequence of the interplay be- 
twee n turbulence, gravity , and feedback in the first place 
fe.g.. lHopkins et al.ll201ll ). 
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always stable against perturbations. However, if the SFR 
changes strongly with metallicity the equilibrium can be- 
come unstable. In this case the perturbation initially grows, 
i.e., gas is enriched or depleted, but eventually Jg- becomes 
small (e.g., when all the gas is molecular) and the system 
approaches a new equilibrium. We note that this derivation 
ignores changes to the gas masses (or stellar masses), and 
hence, it does not prove that Z cq is the metallicity of an equi- 
librium state that is also stable against these perturbations, 
which is what we show next. 

Consider a galaxy that resides in a halo with mass Mhaio 
and that accretes gas and stars according to the Equations 
Q and ((5]). The accretion rates are governed by gravitation 
and not affected by changes in the gas fraction or the metal- 
licity of the galaxy. Also, the scale radius of the gas disk 
(proportional to the virial radius of the halo) is not affected 
by such changes (at least in the context of the model of Sec- 
tion [2J. We assume that the star formation in the galaxy is 
described by a power-law relationship 



APPENDIX A: EQUILIBRIA AND STABILITY 
ANALYSIS 

The time derivatives of the metallicity Z, metal fraction fz, 
gas fraction f g , and stellar fraction f s are 

(A2) 

here X £ {Z, g,s} and we introduced the functions gz, g g , 
and g s . 

Combining Equation (|A1|) with Equations Q and (fl"2|) 
we obtain 



Z = 



— [2/(1 - R)(l - C)M»,for m ~(Z- ZlGM)Mg,in] 



M 



M g , in 



[y(l-R)(l-0r-(Z-Z 1GM )]. 



(A3) 



The second equality assumes M g ,i n 7^ and we introduced 
r = M^form/A/g.in- Clearly, aside from trivial cases such as 
Af»,f rm = and Z — Zigm, the metallicity can only be 
constant if the equilibrium condition (|2UI2ip is satisfied, i.e, 
for Z cq = y(l -R)(l- ()r + Zigm. 

Let us assume that the metallicity is perturbed slightly 
relative to its equilibrium value, i.e., Z — Z cq + S with a 
small perturbation 8. Let us also introduce the shorthands 

A = 1 - R + tout, B — 2/(1 — R)(l - 0, G = £in/g,in/b, 

D = (1 — / g ,in)/b to simplify the notation. The evolution of 
S can be obtained from the linearized version of Equation 



5 = 3 



M s ,in 



B^(Z cq) 



(A4) 



Hence, we find that <5 — > if £>J-§ < 1. In particular, if, 
for a given gas configuration, star formation proceeds in- 
dependently of its metallicity, the equilibrium metallicity is 



M., tam = N\f a f Ha (fz)] % 



(A5) 



Here, N is a normalization constant (which includes every- 
thing that does not change when the gas fraction or metal 
fraction changes) and n is the power-law exponent of the 
star formation - gas relation. The assumption that the H2 
fraction of the gas depends on fz and not on Z deserves a 
comment. Equation ((H) shows that the H2 fraction depends 
primarily on the dust optical depth (Equation 1 11[) . The lat- 
ter is proportional to ZE g oc fz- 

Let fx = fx,eq + $x be a small perturbation around 
the equilibrium state (/z, cq , / g , cq , /s, eq ) (see Equations 1221 
and I25p. The evolution of the perturbation is governed by 
the following system of equations 



fx 



9x{fz, fg,fs 



E 



Y£{Z,g,s) 



^ Y O-f (fz,cq,fg,aq,fs,cq)- 

OJY 



If all eigenvalues of the Jacobian Jx,y = have a 

real value less than zero, then a system in the state 
(fz,eq, /s,eq, /s,eq) is stable against small perturbations. 

In the following we use the additional shorthands E = 
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9 In /. 

1 
T 
1 

T 

Jz.s = 



Jz,z 

Jz,g 



(,fz, cq ) and T = A/ halo /A/ halo . We find that 

A 



(fz,eq — CZlGM 
(fz,eq — CZlOu) 



nE 



f Z , c q B fg , c q A fz , c q 
n fz,eq A 



_/s,eq /g,eq Bf g>cq Af, 



Z,cq 



7l£ 



T 
1 

r L/ S ,e 



/Z,cq 

n 



(C — /s,eq) 
(C - / 9 ,cq) + 1 



J 3 ,s = 



Js,Z 



nE 



fz,cq 

1 n 



2^ /s,cq 

_ JL 



(/ S ,eq - O) 



Replacing (/z, cq , /s,eq, /s,eq) with the expressions (l22j) - 
(|24[) allows us to write Jx,y as a function of the following 
parameters: n, A, B, E, Zigm, and r, the ratio between SFRs 
and gas accretion rate, for the given equilibrium state. In 
particular, C and D drop out, which shows that the stability 
of the equilibrium does not depend on the values of / g ,j n or 

tin. 



Jz.z — — 



nEBr 2 A + (nEZi GM A + B - nEB)r + ZlGM 





1 




~ T 


Jz,s 


= 




1 


Jg,z 


~ T 




1 


Jg,g 


~ T 


Jg,s 


= 


Js,Z 




Js,g 




<Js,s 





(nBA 



(Br + Zi GM ){Ar - 1) 
BA)r -nB + nZ 1Gm A - Z IGM A 



Ar - 



nEAr 



(Br + Z IGM ){Ar - 
1 + (n- l)Ar~ 



Ar - 



27 



rnE 



50 (Br + Zi GM )(Ar - 
27 nr 



1) 



50 Ar 



Clearly, J has a singularity at Ar~l as expected from Equa- 
tion (1221) and it has an eigenvalue -1. The convergence to 
the equilibrium solution occurs on the timescale T if the 
other two eigenvalues have real components that are smaller 
or equal to this eigenvalue. The general expression for the 
eigenvalues of J is rather lengthy and we do not reproduce 
it here. Instead, we compute the real part of the eigenvalues 
numerically using the default parameters for A, B and Z\ G m- 
The parameter B depends on £ and thus on halo mass. We 
adopted £ = which results in the most stringent limit on 
stability. 

We show in Fig. lAll the result of our numerical analysis 
for n — 1 and n — 2. The maximum of the real part of the 




iog 10 r 



Figure Al. Maximum of the real part of the eigenvalues of JT 
as function of the ratio between SFR and gas accretion rate. The 
equilibrium is stable if the maximum of the real part of the eigen- 
values is less than zero (horizontal dot-dashed line). Results for a 
linear (n = 1) star formation - H2 relation are shown by the three 
solid lines. The individual lines correspond to E = 1, E = 0.3, 
and E = 0.1 from top to bottom. Dashed lines show the cor- 
responding results for a quadratic star formation - H2 relation. 
The equilibrium is stable for r < 1/A (vertical dotted line) if 
E < 0.5 — 1 (depending on n), i.e., as long as the H2 fraction is 
not extremely sensitive to changes in the metal fraction. 



eigenvalues is less than zero (and hence the equilibrium is 
stable) for any r < 1/A as long as E is sufficiently small (E < 
0.5 — 1). Hence, as long as the H2 fraction is not extremely 
sensitive to changes in the metal fraction, the equilibrium is 
stable. 
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